Fourier analyses of commensurability oscillations in Fibonacci lateral superlattices 
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Magnetotransport measurements have been performed on Fibonacci lateral superlattices (FLSLs) 
— two-dimensional electron gases subjected to a weak potential modulation arranged in the Fi- 
bonacci sequence, LSLLSLS..., with L/S=t (the golden ratio). Complicated commensurability 
oscillation (CO) is observed, which can be accounted for as a superposition of a series of COs each 
arising from a sinusoidal modulation representing the characteristic length scale of one of the self- 
similar generations in the Fibonacci sequence. Individual CO components can be separated out from 
the magnetoresistance trace by performing a numerical Fourier band-pass filter. From the analysis 
of the amplitude of a single-component CO thus extracted, the magnitude of the corresponding 
Fourier component in the potential modulation can be evaluated. By examining all the Fourier 
contents observed in the magnetoresistance trace, the profile of the modulated potential seen by the 
electrons can be reconstructed with some remaining ambiguity about the interrelation of the phase 
between different components. 
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I. INTRODUCTION 

A GaAs/AlGaAs-based high mobility two-dimensional 
electron gas (2DEG) combined with modern electron- 
beam (EB) lithography allows us to fabricate a device 
in which electrons travel through an artificially designed 
environment over a sufficiently long distance before be- 
ing scattered by impurities fi A prototypical example of 
such devices is a lateral superlattice (LSL), where a 
periodic modulation of electrostatic potential is intro- 
duced to a 2DEG. The periodic modulation can be ei- 
ther one-dimensional (ID) or two-dimensional (2D), ex- 
emplifying the versatility of LSLs compared with more 
conventional vertical superlattices. The planar geom- 
etry of a LSL is well-suited for experimental studies 
of the transport properties, and therefore the magne- 
totransport of both ID and 2D LSLs, with the mag- 
netic field B applied perpendicular to the 2DEG plane, 
have been the subject of extensive studies se e, e.g.. Refs. 
11^3,45,6X1100 IL 12 13 14 15 and Refs. [llll^liiillO 
(and references therein) for ID and 2D LSLs, respec- 
tively] . 

In a ID-LSL, a number of phenomena are known 
in the magnetotransport, including low- field positive 
magnetoresistance (PMR),^ commensurability oscilla- 
tion (CO) originating from the commensurability be- 
tween the cyclotron radius and the modulation period;^ 
quantum interference of closed orbits and geomet- 
ric resonance of open orbits. ^» The periodic modula- 
tion in a ID-LSL can be replaced by a quasiperiodic 
one without appreciable technical difficulties. Quasiperi- 
odic structures are characterized by long-range order 
without translational symmetry, and constitute intrigu- 
ing intermediate states situated between perfect pe- 
riodicity and disorderi^ A wealth of exotic proper- 
ties, such as localized or critical wave functions and a 
Cantor-set spectrum, are theoretically predicted,— i^^i^ 
with varieties of experimental systems devoted for 
their exploration. 2^^i21^iS9^i^^^^^^^^^ It is 



therefore of great interest to perceive how the magne- 
totransport phenomena observed in a periodic ID-LSL 
are altered (or remain unaltered) in a quasiperiodic ID- 
LSL, as well as to search for phenomena peculiar to the 
quasiperiodic systems. The Fibonacci sequence typifies 
quasiperiodic systems in one-dimension. The present au- 
thors have been studying the magnetotransport of LSL 
samples with the potential modulation arranged in the 
Fibonacci sequence — the Fibonacci LSL (FLSL). In 
previous publications, we briefly reported complicated 
CO and geometric resonance of open orbits observed in 
FLSLs)^ as well as the enhancement of zero-field resis- 
tivity followed by strong negative magnetoresistance for 
FLSLs with short-length-scale modulation approaching 
the Fermi wave length4i In the present paper, we fo- 
cus on the quantitative Fourier analysis of CO. We will 
show that the magnitudes of individual Fourier compo- 
nents that constitute the potential modulation can be 
determined by the analysis. Detailed knowledge of the 
potential profile thus obtained forms the basis for un- 
derstanding the electronic properties or phenomena in 
FLSLs to be explored in the future studies. 

The paper is organized as follows. In Sec. [TTl the CO 
for a simple sinusoidal periodic LSL and its extension to 
systems with multiple Fourier components are briefly re- 
viewed, in order to prepare for later application to the 
analysis of FLSLs. Details of the samples and the exper- 
imental setup are described in Sec. lIIII The main subject 
of the paper is presented in Sec. IIVI where we delineate 
the prescription to single out CO corresponding to a par- 
ticular Fourier component from the magenetoresistance 
trace of a FLSL, and to deduce the amplitude of the po- 
tential modulation for the component; we further make 
an attempt to reconstruct the potential profllc by ex- 
ploiting the phases of the Fourier components inferred 
by examining the origin of each component. Discussion 
on the limitation in the accuracy of the potential profile 
thus determined is given in Sec.|Vl followed by conclud- 
ing remarks in Sec. IVII 
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II. COMMENSURABILITY OSCILLATION IN 
PERIODIC LATERAL SUPERLATTICES 

In a periodic ID-LSL, the oscillatory component of the 
resistivity in the CO, 5p%^{B), is well understood quan- 
titatively. We take the cc-direction to be the direction of 
the modulation V{x) and also the direction of the cur- 
rent for the resistivity measurement, and the a;?/-plane 
to be the plane of the 2DEG with B applied along the 
2;-axis, throughout this paper. For a sinusoidal modu- 
lation V(x) — VgCos{gx) with period a=2'n j g and with 
the amplitude Vg much smaller than the Fermi energy 
E-p, 5p^^{B) has been experimentally shown to be accu- 
rately described by a simple analytic formula(^i^ 

'''-^''^ =,A(-^)gV^BM29R.), (1) 
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with ctq = riee^Ttn/m* for (u;crm)^=(/xi3)^»l. This, with 
the Bessel function replaced by its asymptotic expression 
Jo(X)~^2/7rX cos(X - 7r/4), coincides with Eq. (P) ex- 
cept for non-oscillatory term linear in B and the damping 
factors. The scattering damping factor can be incorpo- 
rated within the semiclassical picture by the Boltzmann 
equation approach. The thermal damping factor results 
from standard treatment of the blurring of the edge of the 
Fermi-Dirac distribution function. Virtually the same re- 
sult (with A{T/Ta) but without A{n / ^^B)) is arrived at 
by quantum mechanical calculations that treat V{x) by 
the first-order perturbation theory4^i^ 

By fitting Eq. ([T]) to an experimentally obtained CO 
trace, unknown parameters of the LSL sample under in- 
vestigation can be deduced. It can readily be seen from 
Eq. ([1]) that the oscillation is periodic in with min- 
ima taking place at the conditions. 
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Po=Px.(0), A{X)=Xlsm\i{X), k^Tg^{ll2^){kylg)huj,, 
m* the electron effective mass (0.067me for GaAs), 
uJc=eB /m* the cyclotron angular frequency, Rc=hkp/eB 
the cyclotron radius, fcF=\/27me the Fermi wave number, 
and fi and rie represent the mobility and the density of 
the electrons, respectively. The factor A{Tr/^^B) in Eq. 
P]) is the damping factor arising from the scattering of 
electrons out of their cyclotron orbit, with fi^ an ap- 
propriate mobility usually identifiable with the quantum 
mobility fiq that describes the damping of the Shubnikov- 
de Haas (SdH) oscillation. Apart from the thermal 
and scattering damping factors A{T/Ta) and A{tt/ fi^iB), 
Eq. ([T]) can be deduced from a physically transparent 
semiclassical picture4i Electrons acquire a drift veloc- 
ity Vii^y{x)——E{x)/B in the simultaneous presence of an 
electric field E{x)—{l/e)dV{x)/dx and a magnetic field 
B. The average of the drift velocity over one cycle of the 
cyclotron orbit 
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(n = 1,2,3,...). 
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From the frequency Bg, the wave number g='KeBg/hkF 
of the potential modulation (if unknown) is revealed us- 
ing the Tie determined either by the Hall resistivity or by 
the SdH oscillation. The parameters Vg and p-w are ob- 
tained from the oscillation amplitude. Fitting the exper- 
imental values of \6p^^{B) / poA{T /Tg) \ at the extrema to 
the function C/ s'mh{By^ / B) employing C and Sw as fit- 
ting parameters, we get Vg—y^C/{Byj^g) and fivj—ir/B^, 
where 7 is evaluated by Eq. ([2|) using experimentally ob- 
tained parameters /i and rig. 

Extension of the above argument to general poten- 
tials of the form V{x)—J2g cos{gx + (j)g) is straight- 
forward. The resulting conductivity was shown, using 
the quantum mechanical Kubo formula, to be simply 
additive*^ it can also readily be shown by the semiclassi- 
cal derivation presented above by noting that the cross- 
terms in the square of the cycle-averaged drift velocity 
"^(^^o ) =Z]g /eB) sin(5ra;o + 0g ) Jo (5-Rc ) vanish upon 
averaging over xo. The total oscillatory part of the CO 
thus reads 
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oscillates with B reflecting the commensurability be- 
tween a and i?c, with the amplitude depending on the 
guiding center position xq. Here, Jo{X) represents 
the Bessel function of order 0. The conductivity de- 
riving from the commensurability effect, 5a yy^ is ob- 
tained by further averaging the square Wd,j,^(a;o) over 
X{j (averaging over one modulation period suffices for 
a periodic LSL), and substituting the resultant {vj^'^) 
in the Einstein's relation, 5ayy=e'^{'m* /ttTi'){t^((v^'^)), 
with T^={e/m*)fi the (momentum relaxation) scattering 
time. The corresponding resistivity is obtained by invert- 
ing the conductivity tensor: Spxx/po=i^cTni)'^5ayy/aQ 
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with Sp^x given by Eq. ([T]). The phase (jjg is equivalent to 
the shift of the guiding center by ipg/g, and therefore does 
not appear in Eq. ([5]) obtained by the averaging over xq; 
the information on the phase is lost in the resistivity. For 
a (non-sinusoidal) periodic potential that includes higher 
harmonicSj^ the g in the summation contains only in- 
teger multiples of the fundamental component go, i.e., 
g—Xgo (A=l,2,3,...). However, Eq. ([5]) is apphcable to 
more general potential that involves g's with the ratio A 
given by irrational numbers, as is the case in our present 
FLSL samples. Once experimentally measured (5/9^°* is 
decomposed into its constituent Sp^^^s, the correspond- 
ing Vg's can be obtained following the prescription for a 
sinusoidal modulation described above. This is the main 
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FIG. 1: (Color online) Schematic illustrations of the FLSL 
samples, (a) and (b) represent L-type and S-type samples, 
respectively. 



TABLE I: List of samples. (Aires) denotes the average dis- 
tance between the centers of adjacent resist slabs. 

Sample ID L (nm) 5* (nm) (Aa::rcs) (nm) Type 
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subject of the present paper that will be demonstrated 
in what follows. 



III. EXPERIMENTAL DETAILS 

The Fibonacci sequence, LSLLSLSLLSL... , is gen- 
erated by first preparing two unit lengths, L and S with 
the ratio L/S=t={1 + x/5)/2=l. 61803... (the golden 
ratio), and then repeating the inflation rules, S L 
and L — > LS, starting from a single S. In a FLSL, 
modulation arranged in the Fibonacci sequence is in- 
troduced to a 2DEG. The samples are fabricated from 
conventional GaAs/AlGaAs single heterostructurc 2DEG 
wafers (/i=50-100 mVVs, ne=1.7-3.0x 10^^ m'^, varied 
by LED illumination) with the heterointerface residing 
at the depth of 90 nm from the surface. The electro- 
static potential modulation is introduced by exploiting 
the strain-induced piezoelectric effect,^" the strain be- 
ing generated upon cooling by the slabs of negative EB- 
resisf*^ placed on the surface selectively on either of the 
L-site (L-type FLSL) or S'-site (S'-type FLSL), as de- 
picted in Fig. [TJ Seven i-type and three S'-type samples 
with varying unit lengths, as tabulated in Table IH are 
examined. The accuracy of the actual width of the resist 
slabs is estimated from the scanning electron micrographs 
to be roughly within ten percent of the ideal values. 
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FIG. 2: (Color online) Magnetoresistance traces of a FLSL 
(solid line) and the adjacent plain 2DEG (dotted line) taken 
at 4.2 K from sample 70L (left axis). Magnetoresistance trace 
for a periodic LSL (a=184 nm, Vg=0.3 meV) is also shown 
for comparison (dash-dotted line, with negative offset, right 
axis) . The inset schematically depicts the configuration of the 
sample. 

A Hall-bar pattern with two sets of voltage probes are 
employed for the samples, as illustrated in the inset of 
Fig. [21 allowing for the simultaneous measurement of the 
unpatterned plain 2DEG adjacent to the FLSL for ref- 
erence. The width of the Hall bars and the center-to- 
center distance between the voltage probes are 14 ^m 
and 40 /im, respectively. Thus the FLSLs contain from 
140 (sample HQS') to 290 (sample 40L) segments of the 
S'-site. The main panel of Fig. [2] shows a typical mag- 
netoresistance trace of a FLSL, along with that of the 
adjacent plain 2DEG. The measurement is carried out 
by standard low-frequency ac lock-in technique. We have 
checked that the microfabrication process to introduce 
the modulation does not degrade the mobility and also 
does not alter the electron density (as is evident by com- 
paring the SdH oscillations observed in the two traces for 
B>Q.b T). In the FLSL, the CO is observed for S>0.05 
T, following the low-field PMR. The complicated non- 
monotonic nature of the CO in the FLSL is immediately 
evident by the comparison with the CO of a periodic LSL, 
also plotted in Fig. [2] for reference. Quantitative analysis 
of the CO in FLSLs will be presented in the following 
section. 



IV. FOURIER ANALYSES OF THE 
COMMENSURABILITY OSCILLATION 

A. The separation of Fourier components and the 
analyses of the CO amplitudes 

In this section, we describe, taking sample 70-L 
as an example, the procedure to separate out 
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FIG. 3: (Color online) (a) The plot of {d^ /dB^){Ap^^/po) vs. 
1/B. (b) Harmonic components of (a) obtained by performing 
Fourier band pass filters. Corresponding windows are shown 
in Fig. |4] The left axis is for /s, and other traces are offset 
for clarity. 



Fourier components from a magnetoresistance trace 
^Pxx{B)/ po^{pxx{B)- po)/ po as shown in Fig.[21 and to 
analyze the CO amphtude of each Fourier content thus 
singled out. The first step is to subtract the slowly- 
varying background to obtain the oscillatory part. In 
the case of a periodic LSL, the background can readily 
be defined as the average of the upper and lower enve- 
lope curves;^ since the CO amplitude decreases mono- 
tonically with decreasing magnetic field, the envelope 
curves obtained as spline curves connecting maxima or 
minima, and hence their average, are also monotonic. A 
similar approach is not applicable for FLSLs having non- 
monotonic envelope curves. We, instead, make use of the 
numerical differentiation by B that effectively operates 
as a high pass filter. The second derivative successfully 
eliminates the slowly-varying part, as can be seen in Fig. 
|3l^a) which shows {d? / dB'^)[/S.pxx{B) / po] plotted against 
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FIG. 4: (Color online) The Fourier spectrum of the trace of 
{d^/dB'^){Ap^^/po) vs. 1/B shown in Fig.[3](a). The shaded 
areas indicate the windows employed in the band-pass filters 
in Fig. H 



The Fourier transform of Fig.[3]Ja), presented in Fig. 21 
displays discrete peaks. This indicates that the beating of 
the oscillation in Fig.[31[a) can be viewed as resulting from 
the superposition of a number of components periodic 
in l/B. A peak originating from the SdH effect is also 
observed. The positions of the principal peaks marked 
in Fig. 3] by fj (j=5,4,3,2 and a small trace for are 
mutually related by fj/fj+i = t, namely, the positions 
of adjacent peaks are scaled by the golden ratio r. The 
other two minor peaks bear relations /2.-2 = 2/4 and 
/o,2 = 2/3. These peaks can be explained in terms of the 
reciprocal lattice for a ID Fibonacci sequence. 

It is well known that quasiperiodic structures show a 
set of discrete sharp diffraction spots despite the absence 
of definite periodicity^Si The reciprocal lattice of a ID- 
Fibonacci lattice can readily be found by noting that the 
Fibonacci sequence can be generated also by the pro- 
jection of a 2D-square lattice, with the lattice constant 
yl + r^S, onto a line having the slope equal to 1/r; the 
reciprocal lattice, or equivalently the position of peaks 
in the Fourier transform, is given by gm.n=i'^'^/ S)fm,n 
with 



fm,n — 



(6) 



and m, n integers;^ Although in principle gm,n's fill up 
the entire ID reciprocal space, only a small fraction of 
them with smaller values of |r7i|, \n\ have practically ef- 
fective amplitudes of the structure factor, the detail of 
which depending on the actual profile of the potentiali^ 

The above-mentioned principal sequence coincides 
with Eq. ^ with two successive Fibonacci numbers. 
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i^a-j and n=i^2-j 
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if we translate the frequency Bg into S/2/S.Rc, the inverse 
of the corresponding increment in the cyclotron diame- 
ter 2ARc=2hkF / eBg taking S as the unit of the length 
(the top axis in Fig. 0]). Here the Fibonacci numbers 
F„ are defined by Fq^O, Fi^l, and F„+2=F„+i + F„ 
for both positive and negative integers n, and can be 
exphcitly written as F„ = [r" - (-t)-"]/\/5 (the Bi- 
net formula— 1^). In Eq. ([7]), we made use of an identity, 
T^=TFn + Fn-i. The two minor peaks are also expressed 
in terms of Eq. 

The peak at fm,n corresponds to the modulation wave 
number g—{2-K/S)fm.n or the modulation period S/ f„i,n- 
Therefore Fig.|4]reveals that the profile of the modulation 
is mainly composed of the superposition of multiple in- 
commensurate components having the periods S/ fj , the 
periods successively scaled by an irrational number t. 
Physical interpretation of these periods will be presented 
in the following subsection. 

The trace in Fig. [3] (a) can be decomposed into its 
Fourier components by performing a numerical Fourier 
band-pass filter. The resultant components obtained by 
applying the windows that encompass the corresponding 
peaks (shaded areas in Fig. [5]) are displayed in Fig.[3](b). 
The next step is to restore the resistivity Sp^^ from its 
second derivative. This can be readily done by numeri- 
cally integrating the second derivative (the replot of Fig. 
[3] (b) against B) by B twice. The decomposed Fourier 
components of the CO (and also the SdH oscillation) thus 
obtained are plotted in Fig. \5\ The summation of all the 
observed components (thick dashed curve) restores the 
oscillatory part of the original magnetoresistance trace 
(thick solid curve) with high fidelity, indicating that we 
are on the right track. As shown in Fig. [6] (a), the plot of 
the index of the extrema in each CO component versus 
their position in 1/B falls on a line with the slope Bg that 
coincides with the peak position of the FFT spectrum 
(the lower axis in Fig. 2]) within experimental error, and 
with the ordinate intercept '~-^l/4, as expected from Eq. 
dU, thereby confirming that adequate CO components 
are extracted. Here we assigned half-integer indices to 
the maxima. 

The final step is to deduce the modulation amplitude 
Vg from the amplitude of CO following the procedure de- 
scribed in Sec. [ID The values of \SpS^^{B)/ poA{T/Tg)\ at 
the extrema are plotted in Fig. [H] (b) for the three major 
components, /2, /a, and f^. As can be seen, the function 
C/ smh{Byj/B) describes the behavior of the normalized 
amplitudes quite well. Fitting by the two parameters C, 
Sw can be carried out without difficulty for the compo- 
nents with large enough amplitudes, such as /a and /4 in 
the present case. The scattering parameter p,^—T:/B^ is 
expected not to vary among components. This is found 
to be the case for /a and fi] the values obtained by the 
fitting, /iw=10.2 and 11.3 m^/Vs, respectively, roughly 
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FIG. 5; (Color online) Experimental magnetoresistance trace 
(thick solid line), the Fourier components extracted by the 
procedure described in the text (thin solid lines, with down- 
ward offsets), and the addition of the Fourier components 
(thick dashed line, with upward offset). 



coincide (within experimental error ± ~0.5 m^/Vs). For 
components having smaller amplitudes as J2 in Fig. [7] or 
other components not plotted in the figure, the two pa- 
rameter fitting is found to be difficult to perform reliably. 
In such cases, we resort to the fitting by one parameter 
C with a fixed B^ taken from the component with the 
largest amplitude (/a in the present case), relying on the 
invariability of the /iw Since only the ratio CjB^ is re- 
quired for determining the Yg , possible small error in B^ 
does not affect the resulting Vg very much, being com- 
pensated for by the responding change in C obtained by 
the fitting. (Note that C/ s\tl\\{B^ / B)^C B / B^ for smaU 
B^/B.) 

The values of deduced Vg are plotted against 
f—{S/2'K)g, namely, the values of fj or fm,n of each com- 
ponent given by Eqs. ([7]) and ©, in Fig. [7] (a) for all the 
FLSL samples examined in the present work. As can 
be seen, dominant components are /4, /a, or /2 in most 
samples. Minor components /2.-2 and /o,2 are occasion- 
ally detected for samples with longer unit lengths. An 
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FIG. 6: (Color online) (a) Plots of the index versus position 
(in the inverse magnetic field) of the extrema of CO Fourier 
components shown in Fig. [5] The lines are linear fit to the 
data points, (b) Selected plots of the normalized amplitude 
versus position of the CO extrema. The curves are the fit 
to C/sinh(_Bw/B). Maxima and minima are plotted by open 
and solid symbols, respectively. 

important factor that controls the modulation amplitude 
Vg of the component g is the corresponding length (the 
period) 2-^ j g—Sj ] relative to the depth d of the 2DEG 
plane, owing to the exponential decay oc exp(— gd) along 
the depth of the effect exerted on the surface. In Fig. 
[7](b), we replot Vg against Sjj. Dominant components 
for each sample are found to roughly follow the function 
that fits the modulation amplitude of periodic LSLs with 
period a—i-K/g at ne=2.8xl0^^ m^^, 

Vg = (-e)0ie'sVeTF™ (8) 

with 01=43.8 meV, d^lOO nm (the effective depth of 
the 2DEG that takes account of the average distance 
of the 2DEG wave function from the heterointerface), 
and etf™ representing the dielectric function of 2DEG 
in the Thomas-Fermi approximation using the Fang- 
Howard wave function.— 



FIG. 7: (Color online) (a) Amplitude of the potential mod- 
ulation Vg deduced from the CO amplitude for Fourier com- 
ponents specified by /. (b) Replot of (a) against the length 
S/ f corresponding to the Fourier component /. Dotted curve 
shows amplitude for periodic LSL, Eq. Solid and open 
symbols are used for L-type and S-type samples, respectively. 
Lines are guide for eyes and connect only the data points be- 
longing to the principal sequence, fj. 



Now that we have determined Vg, the profile of the 
modulated potential can be reconstructed with further 
supply of the knowledge of the phases, (j)g. Unfortu- 
nately, however, the information of the phases does not 
manifest itself in the CO, as mentioned earlier. We can 
only surmise the phase of Fourier components by looking 
into their origin, which is the subject of the subsequent 
subsection. Note however that in principle any phase 
relation between an arbitrary pair of components with 
irrational wave-number ratio can be found at some place 
in the sequence, provided that the sequence is infinitely 
long. 
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Edges 




FIG. 8: (Color online) Center: A diagram that sketches the Fibonacci sequence (L and S are presented by dark (blue) and 
bright (red) rectangles, respectively) with the vertical lines above and below indicating the locations of characteristic features: 
the centers of the segments not discriminating their types {{L,S}), the centers of the L segments (L), the centers of the S 
segments (5*), the centers of the L areas counting the two consecutive L's, LL, as a single area {{L,LL}), the centers of 
LL's only (LL), the centers of single isolated L's surrounded by S's on both sides {SLS). The boundaries between L and S 
segments, denoted henceforth as |S| with the |'s symbolizing the edges, are displayed in the top panel (Edges). The edges \S\ 
can be categorized into two groups, namely the edges between LL and S {\LL\) and between isolated L and S {S\L\S), which 
are shown by solid and dashed lines, respectively. Right: Schematic representation of the Fourier component (the "average" 
periodic modulation) originating from the characteristic feature depicted to the left. Left: The inflation and the deflation of 
the Fibonacci sequence. 



B. The phase of the Fourier components 

In this section, we make an attempt to interpret the 
origin of the Fourier components observed in the CO. An 
important property of the Fibonacci sequence to be con- 
sidered in this respect is its self-similarity. As depicted in 
the left panel of Fig.[8l the operation of the inflation rule 
(dividing an L segment into two segments with the ratio 
of t:1 to get new L and S, and renaming an S segment 
as L) transforms the Fibonacci sequence to the Fibonacci 
sequence itself but scaled down by the factor l/r, while 
the deflation rule (the inverse process of the inflation) 
yields the Fibonacci sequence scaled up by the factor t. 
The self-similarity implies the presence in a FLSL of mul- 
tiple length scales scaled to each other successively by the 
factor r, where we refer by the term "length scale" to the 
typical distance between characteristic features that take 
place repeatedly. The ratio /j+i//j=r between the ad- 
jacent major components fj can thus be interpreted in 
terms of the self-similarity. 

An obvious example of such repeatedly occurring fea- 
tures is simply the segment S. Recalling that the CO 
results from the averaging of the square of the drift 
velocity over the position of the guiding center xq, it is 
the spatially averaged distance between the features that 



is reflected to the frequency Bg of the CO. As depicted 
in Fig. [51 the distance between the centers of the ad- 
jacent S segments is either 2L + S=t^S oi L + S=t^S, 
reflecting the fact that L segments appear in either of the 
two consecutive {LL) or single isolated {L) form. Not- 
ing that the former (the latter) takes place in one-to-one 
correspondence with the L (S) of the Fibonacci sequence 
retreating one generation by the deflation process, their 
relative abundance approaches t:1 for a long enough se- 
quence. Therefore, the average distance between the cen- 
ters reads, making use of the relation t^=t + 1, 

accounting for the component in the CO, which is 
one of the dominant components in many of the FLSL 
samples (see Fig. EJa)). Since LL is composed of seam- 
less EB-resist in L-type samples (seamless open gap in 
S'-type samples), S/f^ represents as well the average dis- 
tance between adjacent slabs of the resist or between ad- 
jacent free spaces. The values of (Axres)='S'//3 for the 
samples examined in the present study are tabulated in 
Table m Other principal components, fj, can roughly 
be viewed as deriving from the inflation or deflation of 
/a. More specifically, /i, /2, /4, and /s are attributable 
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TABLE II: The a;-coordinate of the fc-th appearance of char- 
acteristic features defined in Fig. [H] 



location (in unit S) 




K/ ji — yLj Zt J\i + rrac[[fi; + ^) n\ ~r rrac(^/c/ rj| 


Xk 


fc//2 - r/2 - Frac(fc/r) 


4 


fc//3 - 1/2 - rFrac(fc/r) 


Xk 


k/fi - r^Frac(fc/r) 


Xk 


fc//5 + r/2 - Frac(fc/r) 


4 


fc+//o,2 - rFrac(fc+/2r) - 0+ 


\LL\ 

4 


fc+//2,-2 - Frac(fc+/2T) + (-l)'^r 


S\L\S 
Xk 


fc+//_2,4 - r^* Frac(fc+/2T) + 0+r 



TABLE III: The spatially averaged location for of x^s given 
in Table [n] for infinitely long sequence. 



location (in unit 5) 



{4) 
(4) 

i4') 
(4"'> 



k/fi - 1/r 
k/h - 
k/h - r'/2 
k/U - tV2 
k/fs - rV2 
(fc + l/2)/2/3 - rV2 
(fc + l/2)/2/4 - rV2 
(fc + l/2)/2/5 - rV2 



to the centers of both types of the segments, the L seg- 
ments, the LL's, and the isolated L's, respectively (see 
Fig. [5]). The minor components fm,n not categorized in 
fj can be ascribed to the edges of the resist slabs, which 
basically bear the frequency twice of the corresponding 
major components. 

To gain insight into the phase of the modulation, we 
need to take a closer look at the position of the charac- 
teristic features. Since the fc-th appearance (fc=l,2,3,...) 
of the L segment {S segment) takes place at the a^- 
th (&fc-th) place in the whole sequence, where ak=k + 
Int(fc/T) and &fc=2fc-|-Int(fc/r) represent a pair of Beatty 
sequences^ with Int(x) the integer part oix, the position 
of the fc-th center of the S segment, indicated by (red) 
vertical lines in the central panel of Fig. [51 is given by 



xf = (bk - k)TS + kS-^ 



k fk 

— rFrac — 

h \T 



(10) 

with Frac(x)=x — Int(x) denoting the fractional part of 
X. The positions of other features can be written down 
along the same line, exploiting the correspondence of 
the feature to L or in the inflated or deflated Fi- 
bonacci sequence, and tabulated in Table |TT1 employ- 
ing S as the unit of the length. We used the notation 
fc+ = fc + [1 — (— 1)'^]/2 in the description of the minor 
components, which takes the same number for an odd fc 
and the subsequent even fc. 



Since t is an irrational number, Frac(fc/r) and 
Frac(fc"''/2r) take the value uniformly distributed in the 
interval [0, 1), fluctuating with fc around the central value 
1/2. The distance between the features under consider- 
ation, seen by the electrons traveling in the cyclotron 
orbit, fluctuates accordingly with the guiding center po- 
sition xq. In the context of the CO, we only need the 
average of the value over xq. Therefore the fractional 
parts in Table HIl can be replaced by their average 1/2, 
resulting in the spatially averaged locations (x) tabu- 
lated in Table [1111 For the components deriving from the 
edges, we noted that (—1)'" vanishes on averaging and 
hence fc+^fc-|- 1/2, and also used the relations /o, 2=2/3, 
/2 -2=2/4, and /_2,4=2/5. 

For the reconstruction of the potential profile, it is nec- 
essary to further relate the characteristic features dis- 
cussed thus far to the potential modulation they create. 
This requires detailed knowledge of the strain introduced 
into the GaAs/AlGaAs wafer by the resist slabs, which 
is beyond the reach of the present study. Here we sim- 
ply assume that the potential energy becomes higher at 
the 2DEG areas directly beneath the resist slabs, namely 
at i-sitcs (5-sites) for i-type (S'-type) samples, which 
are probably under compressive strain. Accordingly, 
the Fourier components including their phases are writ- 
ten, referring to Table [ml as ±Vj cos[27r/j(x/S' -I- tV2)] 
for j=2, 4, 5, and TV3 cos[27r/3(x/S' + t"^ /2)] for j=3, 
with the upper (lower) sign for L-type (S'-type) sam- 
ples. For the component /i which includes L-sites and 
5-sites alike, we presume that the L-sites dominate the 
sign: ±Vi cos[27r/i(x/S' + l/r)]. For edges, we as- 
sign the sign opposite to the area in between, namely, 
±Vo,2 cos[47r/3(x/S' -I- r/4)] and TV2-2 cos[47r/4(x/S' - 
1/4t)] (the component /-2,4 was not observed in the 
present experiment). 

The modulation profile is obtained by adding up all 
the components observed in the experiment. Examples 
of the profiles thus reconstructed, as well as the con- 
stituent Fourier components, are displayed in Figs. [51 and 
Unifor L-type (50L and 70L) and S'-type (90S and llOS) 
samples, respectively, using the phases assumed above. 
It is intriguing to note that the superposition of a small 
number of incommensurate periodic sinusoidal modula- 
tions basically reproduce the quasiperiodic profile. Be- 
cause of the low-pass filtering action of the attenuation 
exp(— f/(i)=exp[— (27r//S)(i] mentioned earlier, higher fre- 
quency components (components fj with smaller j), i.e., 
the components with shorter length scale, are damped 
heavier for samples with smaller S. This is reflected in 
the difference of the profile between samples, resulting in 
rather dull line shape for samples with smaller S when 
plotted against x scaled by the unit length S (see the 
top axes). The profile of S-type samples is qualitatively 
similar to that of L-type samples but laid up-side-down. 
Quantitatively, however, the modulation amplitude of S- 
type samples is much smaller than that of corresponding 
L-type samples with the same set of length units L and 
S (compare samples 90L and 90S in Fig. [7]). This is 
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FIG. 9: Potential profile of FLSLs (solid lines) reconstructed 
by the superposition of the Fourier components (dotted lines) 
with the amplitudes determined by the analysis of the CO 
and the phase obtained from the consideration of the origin 
of the components, for L-type samples SOL (a) and 70L (b). 
Note the difference in both vertical and horizontal scales. The 
horizontal axis is scaled so as to keep x/S (top axis) the same 
for both figures. The shades indicate the area in the 2DEG 
plane located right beneath the resist, namely the L-sites. 



attributable to 1/r times smaller coverage by the resist 
slabs of the surface that exerts much smaller strain to 
the underlying GaAs/AlGaAs 2DEG wafer. 



V. DISCUSSION 

What we have shown in the previous sections demon- 
strates that the Fourier analysis of the CO serves as 
a useful tool to probe the complicated profile of the 
modulated potential introduced into a 2DEG. FLSLs 
contain a number of Fourier components and therefore 
provide archetypal proving ground for the technique. 



FIG. 10: Similar to Fig. |9] for S-type samples 905 (a) and 
llOS (b). The shades indicate the S'-sites. 



The method, however, can be applied in wider range 
of systems that are comprised of multiple Fourier con- 
tents, e.g., periodic systems that retain higher harmon- 
ics. The technique is, in a sense, analogous to diffrac- 
tion techniques that determine the crystal structures. 
In fact, Eq. ([5]) with Eq. ^ closely resembles the Pat- 
terson function^'' P{x)—'E,g{\Vg\'^/2)cos{gx) by replacing 
x—2Rl—2Rc — Tr/2g, and therefore can roughly be viewed 
as the autocorrelation function of the modulated poten- 
tial measured by the cyclotron diameter with offset, 2R*. 
In this section, we discuss the limitations of the CO as a 
technique to determine the "structure" of the modulated 
potential. 

Firstly, the information on the phase of the Fourier 
components is lost from the CO, as has been discussed so 
far. Again this resembles the situation in the diffraction 
techniques; the phases are unavailable in the diffraction 
pattern, which records only the intensity of the diffrac- 
tion spots. In both cases, the phases have to be estimated 
from the symmetry of the system or from the physical 
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origin of the component, as has been done in the preced- 
ing section. This can be a major source of uncertainty 
in the obtained reconstructed structures. Note, however, 
that the knowledge of the intensity of Fourier contents 
(without their phase interrelation) suffices for their ex- 
planation in many phenomena, a representative example 
being, of course, the CO itself. 

Secondly, because of the thermal damping factor 
A{TlTg) in Eq. ^ that decreases with increasing 5, the 
CO is less sensitive to components with larger g, i.e., to 
shorter length scales. Therefore higher frequency com- 
ponents can be overlooked in CO traces, especially in 
those taken at higher temperatures. In principle, com- 
plementary use of the geometric resonance of the open 
orbit,— which is a sensitive probe to higher frequency 
components)^ can fill the gap. However, the amplitude 
of the open orbit resonance still remains unexplained 
quantitatively, and its application to the determination 
of the modulation amplitude is yet to be explored. Con- 
versely, the requirement for the electrons to travel at least 
one cycle of the cyclotron orbit before being scattered sets 
an ultimate upper limit, roughly equal to the mean free 
path, in the length scale that can be examined by the CO. 
Practical upper limit will be smaller than this, and will 
be reached when the oscillation resulting from the length 
scale becomes indistinguishable from the slowly-varying 
background (when the peak becomes indiscernible in the 
Fourier spectrum as shown in Fig. [4] due to the low-pass- 
filtering action of the operation of cP/dB'^). Therefore, 
the modulation profile may possibly contain higher or 
lower frequency components that elude the detection by 
the CO. 

Before concluding the paper, we briefly discuss the 
prospective use of a FLSL as a system to experimen- 
tally test the exotic properties in quasiperiodic systems 
that are theoretically predicted i^^'^'^i^'^ The capability to 
experimentally determine the profile of the potential ex- 
perienced by the electrons is obviously of great advantage 
to this end. However, further improvement of the sam- 
ple that enables to control the electron density is desir- 
able, since it will be necessary for many purposes to vary 
the Fermi energy or the Fermi wave number in detail. 
In 2DEG-based devices, this can generally be achieved 
by implementing a uniform gate. Note that in LSLs we 
have to resort to a back gate, since the front surface of 
the wafer is occupied by the grating that introduces the 
modulation. For back-gated LSLs, it was found that the 
modulation amplitude shows strong increase accompany- 
ing the decrease in the electron density by the application 



of the negative gate bias,^^ which might complicate the 
interpretation of the data obtained by the experiments. 
A project to develop LSLs with conducting resist gratings 
that can be biased to compensate for the change in the 
modulation amplitude induced by the back-gate bias is 
now under progress. Another caveat to be born in mind 
is that our FLSLs are basically a two-dimensional system 
where electrons can travel both parallel and perpendicu- 
lar to the one- dimensional modulation. In many theories, 
by contrast, strictly one-dimensional systems are consid- 
ered. It may therefore be necessary for some purposes to 
start the device fabrication from a quantum wire instead 
of a plain 2DEG. 



VI. CONCLUSIONS 

We have shown that the complicated CO observed in 
a FLSL can be decomposed into the constituent Fourier 
components by applying a numerical band-pass filter, 
with the aid of numerical operation of (P /dB"^ for back- 
ground elimination (and later integration by B twice for 
recovery). From the amplitude of the component thus 
isolated, the magnitude Vg of the potential modulation 
responsible for the component can be determined. The 
analysis reveals that the potential profile seen by the 
electrons in a FLSL is basically described by the super- 
position of a small number of incommensurate periodic 
modulations. The periods of the major components are 
given by S/ fj with fj in Eq. ([7]), and are explicable in 
terms of inflation and deflation rules that lead to the self- 
similarity. By examining the physical origin of the com- 
ponents, we can estimate the phase of each component, 
which allows us to reconstruct the modulation profile. 
With the detailed knowledge of the potential profile, we 
believe that a FLSL provides an excellent arena to exper- 
imentally explore exotic phenomena (localization, critical 
wave functions, a Cantor set spectrum, etc.) envisaged 
in the quasiperiodic systems. 
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